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q ■ ABSTRACT 

■ We use two cosmological simulations of structure formation to study the condi- 

tions under which dark matter haloes emerge from the linear density field. Our analysis 
focuses on matching sites of halo collapse to local density maxima, or "peaks" , in the 
initial conditions of the simulations and provides a crucial test of the central ansatz 
of the peaks formalism. By identifying peaks on a variety of smoothed, linearly ex- 
trapolated density fields we demonstrate that as many as ~ 70% of well-resolved dark 
matter haloes form preferentially near peaks whose characteristic masses are similar 
to that of the halo, with more massive haloes showing a stronger tendency to reside 
near peaks initially. We identify a small but significant fraction of haloes that appear 
CNl ' to evolve from peaks of substantially lower mass than that of the halo itself. We refer 

to these as "peakless haloes" for convenience. By contrasting directly the properties of 
these objects with the bulk of the proto-halo population we find two clear differences: 
1) their initial shapes are significantly flatter and more elongated than the predom- 
inantly triaxial majority, and 2) they are, on average, more strongly compressed by 
tidal forces associated with their surrounding large scale structure. Using the two- 
point correlation function we show that peakless haloes tend to emerge from highly 
clustered regions of the initial density field implying that, at fixed mass, the accre- 
tion geometry and mass accretion histories of haloes in highly clustered environments 
differ significantly from those in the field. This may have important implications for 
understanding the origin of the halo assembly bias, of galaxy properties in dense en- 
vironments and how environment affects the morphological transformation of galaxies 
near groups and rich galaxy clusters. 
2 , 
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1 INTRODUCTION 

Within the standard cold-dark-matter cosmology (ACDM) 
non-linear clumps of dark matter (referred to as haloes for 
short) form hierarchically through the merger and accre- 
tion of previous generations of virialized objects. It is within 
these haloes that baryons cool and condense to form galax- 
ies, the visible tracers of the underlying matter distribution 
(White & Rees, 1978). Understanding galaxy formation thus 
requires detailed theoretical predictions for the spatial distri- 
bution of dark matter haloes, their abundance as a function 
of mass, and how these quantities evolve temporally. 
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Attempts to tackle these issues have come from a num- 
ber of directions, the most successful of which has been to re- 
sort to direct N-body simulation of non-linear gravitational 
clustering starting from initial conditions where linear the- 
ory is applicable. Although N-body simulations are power- 
ful probes of non-linear structure formation they are noto- 
riously difficult to implement and provide little theoretical 
insight into the structure formation process. Identifying the 
precursors of dark matter haloes in the linear density field, 
for example, remains an unsolved problem. 

As an alternative to N-body simulations, dynami- 
cal models for the collapse of overdense regions, such as 
the spherical (Gunn & Gott, 1972) and ellipsoidal (e.g., 
Zel'Dovich, 1970; Icke, 1973; White & Silk, 1979; Hoffman, 
1986; Lemson, 1993; Bertschinger & Jain, 1994; Bond & My- 
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ers, 1996) collapse models, can be used to physically mo- 
tivate the conditions under which haloes emerge from the 
linear density field. In the so-called peaks formalism, one 
assumes that gravitationally bound objects form at local 
maxima (peaks) of the linear density field (averaged over 
the appropriate mass scale) which lay above some density 
threshold. It is difficult to trace back the exact origin of 
this model: some of its features were already discussed in 
Doroshkevich (1970) while the whole idea is briefly sketched 
in Peebles (1980) (§26, p. 124). 

The formalism gained great popularity after Kaiser 
(1984) showed that the statistical properties of overdense 
regions can naturally explain the high-clustering amplitude 
of Abell clusters. This was based on the assumption that rich 
galaxy clusters form from rare high density peaks in the lin- 
ear density field, which are much more strongly clustered 
than the underlying mass distribution (see also, Politzer & 
Wise, 1984; Otto et al., 1986; Cline et al., 1987; Peacock 
et al., 1987; Coles, 1989; Mo et al., 1997; Matsubara, 1999) 

The mathematical development of peaks theory was 
laid out in considerable detail by Peacock & Heavens (1985) 
and Bardeen et al. (1986) (hereafter, BBKS). These authors 
studied the properties of local maxima in Gaussian random 
fields, providing a number of useful statistics, such as the 
number density of peaks of a given height, and their clus- 
tering properties. Extensions to the work of BBKS to non- 
Gaussian density fields have also been considered (Grinstein 
& Wise, 1986; Matarrese et al., 1986; Catelan et al., 1988). 

Blumenthal et al. (1984) used the peak model to sketch 
how galaxy formation takes place in a universe dominated 
by cold dark matter. The number density and internal struc- 
ture of haloes forming from linear density peaks was eval- 
uated by Hoffman & Shaham (1985), while their acquisi- 
tion of angular momentum during collapse was discussed 
in a series of papers by Hoffman (1986), Heavens & Pea- 
cock (1988), Steinmetz & Bartelmann (1995) and Catelan 
& Theuns (1996). Calculations of the halo mass function, 
accounting for the "peak-in-peak" problem, have been pre- 
sented by Manrique & Salvador-Sole (1995) while Manriquc 
& Salvador-Sole (1996, 1998) discussed the implications for 
halo mass accretion histories. 

Recently, new attention has been given to the peaks 
formalism in attempts to model the large-scale clustering 
properties of haloes in both real (Desjacques, 2008; Des- 
jacques et al., 2010) and in redshift space (Desjacques & 
Sheth, 2010). The peak model has also been used to investi- 
gate why the clustering properties of haloes of a given mass 
depend on their formation time (Dalai et al., 2008), the so- 
called assembly bias. 

Over the past couple of decades, the body of assump- 
tions that fashion the peaks formalism have become clear, 
and their validity is rarely questioned. The peaks formalism 
appears to work in a statistical sense (e.g., White et al., 1987; 
Weinberg & Gunn, 1990; Park, 1991; Park & Gott, 1991), 
yet the question of whether the peak model for galaxy for- 
mation actually reproduces the sites of halo collapse is still 
largely unverified. Frenk et al. (1988) were the first to ad- 
dress this issue. Using 32 3 -particle N-body simulations of 
relatively small volumes (box sizes of ~ 14 h^ 1 Mpc) these 
authors concluded that dark matter haloes form preferen- 
tially around high peaks identified in the simulations initial 
density field. Later, Katz et al. (1993) challenged these con- 



clusions using a simulation (144 particles) of considerably 
larger volume (~ 50 h^ 1 Mpc box), by noting a rather poor 
correspondence between peak tracers identified in the initial 
linear density field and the haloes that subsequently form. 
A systematic study of proto-haloes (i.e. the regions in the 
linear density field that collapse to form haloes at a given 
time) in a high-resolution N-body simulation (256 3 parti- 
cles in a ~ (85 h^ 1 Mpc) 3 box) was presented by Porciani 
et al. (2002b). These authors concluded that nearly half of 
galactic sized proto-haloes do not contain a linear density 
peak within their local Lagrangian volume. If confirmed, 
these results highlight a general failure of the peak model 
ansatz since they suggest that a simple one-to-one mapping 
between haloes and peaks in the smoothed density field does 
not exist. 

In this paper we revisit these issues using two high- 
resolution cosmological N-body simulations (1024 3 parti- 
cles). We identify the locations of density peaks in our simu- 
lation initial conditions and contrast these directly with the 
actual sites of halo formation, improving upon the short- 
comings of previous work. The remainder of the paper is 
organized as follows. Following a brief description of our sim- 
ulations in §2, we describe our main analysis techniques in 
§3; §3.1 provides a definition of halo mass, and an algorithm 
for locating peaks and matching them to evolved haloes is 
given in §3.2 and §3.3. In §4 we describe our main results, 
which are summarized in §5. 



2 THE SIMULATIONS 

We study the growth of cosmic structure using two fully cos- 
mological N-body simulations of large scale structure. Both 
runs start from Gaussian initial conditions and assume a 
fiat ACDM cosmological model with numerical parameters 
consistent with the WMAP 5-year data release (Komatsu 
et al., 2009). These are: fi M = 0.279, fi A = l-fl M = 0.721, 
o-g = 0.817, n s — 0.96, and a Hubble constant Ho = H(z = 
0) = 73 km s' 1 Mpc. 

Each run follows the evolution of the dark mat- 
ter distribution using 1024 3 particles in a periodic box 
with fixed side-lengths equal to l box = 150 h^ 1 Mpc and 
1200 ftT 1 Mpc. Initial conditions are generated at z — 70 
(150 h,- 1 Mpc box) and z = 50 (1200 h' 1 Mpc box) using 
the Zel'dovich approximation. Integrations are performed 
in comoving coordinates using a lean version of the tree- 
PM code Gadget-2 (Springel, 2005). For this choice of box 
size, resolution and cosmological parameters the particle 
masses are 2.433x 10 s h^ 1 M and 1. 246X10 11 h' 1 M in 
the 150 ft -1 Mpc and 1200 fr 1 Mpc boxes respectively. By 
combining the results from the two runs we are able to 
robustly probe halo abundances for objects with masses 
7.79 x 10 9 /i _1 Mq, as well as gauge the sensitivity of our 
results to numerical resolution. Further details regarding the 
simulations, including a thorough discussion of the simula- 
tion initial conditions, can be found in Pillepich et al. (2010). 

3 ANALYSIS 

In this section we provide a brief discussion of our main 
analysis techniques. This includes a definition of halo mass; 
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a method for generating smoothed density fields, S s ; a choice 
of filter function, and an algorithm for locating peaks in S s . 

3.1 Definition of halo mass 

Dark matter haloes are identified at z = using a standard 
friends-of-friends (FOF) algorithm (Davis et al., 1985). Our 
halo finder has one free parameter, £, which implicitly defines 
the linking length in terms of the mean particle density, n: 
InT 1 / 3 . For well resolved objects, the FOF method selects 
regions enclosed by an isodensity contour oc l/i 3 . In what 
follows we shall adopt I = 0.2, the conventional choice for 
Q = 1 cosmologies. For each halo identified at z — we 
compute the FOF mass, defined as the total mass associated 
with the FOF group, and retain the identity of the most 
bound particle. When necessary, and at z = 0, we identify 
the position of the most bound particle with the halo center. 

One advantage of the FOF method over other halo iden- 
tification techniques is that it imposes no shape restriction 
on the boundary of the system; particles that lie within one 
linking length from any other particle are grouped together 
into a single system, regardless of its final dynamical state. 
A drawback, however, is the tendency of the FOF method 
to occasionally link separate objects in transient dynami- 
cal states, such as fly-by's or ongoing mergers (e.g. Gelb & 
Bertschinger, 1994; Governato et al., 1997), which can lead 
to the appearance of massive haloes through the temporary 
association of two or more low-mass objects. 

We assess the importance of this effect by comput- 
ing, for each FOF halo, the distance between its true cen- 
ter, r m b (defined by the location of its most bound parti- 
cle) and the geometric center of all FOF particles: d g = 
|r m b — i"cm|/i~2oo- As discussed by Bett et al. (2007) and 
D'Onghia & Navarro (2007), the "center-of-mass offset" is 
a simple test of the prevalence of unrelaxed substructure in 
the halo, and can be used to cull the halo catalogs of objects 
far from equilibrium. We have verified that the results pre- 
sented in §4.1 are not particularly sensitive to this selection 
criterion by repeating the analysis for haloes with d fl < 0.1, 
and doff < 0.3. Because of this we shall consider all haloes, 
regardless of the details of their equilibrium state, in the 
analysis that follows. 

3.2 Identification of linear density "peaks" 

For each of our simulations we compute the linear overden- 
sity field <5(x) using cloud-in-cell interpolation; densities are 
computed, in each case, on a 1024 3 grid. The linear overden- 
sity field can then be smoothed on some spatial scale R/ by 
means of convolution with an appropriate window function, 
W(x, Rf). In what follows we mainly restrict our discussion 
to results obtained with a real-space tophat filter, but have 
verified that our conclusions are not particularly sensitive to 
this choice by repeating parts of the analysis with a Gaus- 
sian kernel. 

For small perturbations in the smoothed field, S s = 
5(x, R/), the tophat filter contains a characteristic mass 
M/ ~ (4-7r/3)Rjp within R/ . For each of our runs we com- 

Strictly speaking, the filter radius and mass are related by 
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Figure 1. Comoving number density of peaks identified in our 
simulation initial conditions plotted as a function of filter mass. 
Peaks are identified on a sequence of linearly extrapolated density 
fields smoothed with a spherical filter containing mass TsAf (see 
text for details). Curves extend from Mf = 1.16 X 10 16 /i- 1 M Q , 
the maximum smoothing scale considered, down to Mj = M par t, 
corresponding to the unsmoothed density field. Points are over- 
plotted to highlight the mass range My > 32 X M pa rt, the min- 
imum halo mass identified by our group finder. Squares and cir- 
cles differentiate results obtained with a Gaussian and tophat fil- 
ter, respectively. The tophat filter contains a characteristic mass 
M f = (47r/3)R3p; the Gaussian filter has M f = (2?r) 3 / 2 R3p. Red 
curves (with open symbols) correspond to our 1200 h,- 1 Mpcbox; 
results from our 150 h" 1 Mpc box are shown in blue. The theo- 
retical curves are the peak number densities expected from peaks 
theory (eq. 1) for several different filter choices (see text for de- 
tails). 

pute 5s in equally spaced steps in logM/ of fixed logarith- 
mic width AlogM/ = 0.129. We choose a minimum filtering 
mass equal to the particle mass and perform the convolution 
up to a maximum mass of M/ = 1.16 x 10 16 h^ 1 Mq. This 
results in 60 steps in M/ for our 150 Mpc box; for our 
1200 hT 1 Mpc box, which has a higher particle mass, there 
are a total of 40 steps in Mf. 

Given the smoothed density fields, S s , we next identify 
local maxima, or "peaks", by locating grid cells that are 
denser than all 26 neighboring points. Peaks are identified 
on all smoothed fields and their evolution followed by tag- 
ging the particle at x p nearest to each peak grid point. For 
each peak particle we calculate the associated overdensity, 
5 P — <5s(x p ,R/), the peak height, v v = <5 p /<r(Rf), and the 
properties of the halo (if any) to which the particle belongs. 
In the following we consider all peaks in 8 a and avoid im- 
posing an arbitrary minimum height for their selection. 

Mf oc (1+<5)R 3 , so that for linear perturbations in 5 (i.e., 6 << 1) 
we have Mj- oc R^ 
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With this information we plot in Figure 1 the co- 
moving number density of peaks as a function of smooth- 
ing scale, n p (Mf), for two different choices of filter func- 
tion. Open and closed symbols distinguish results obtained 
from our 1200 /i" 1 Mpc and 150 ft -1 Mpc boxes, respec- 
tively. For the tophat filter points span the mass range 
32 x Mpart < Mf < lO 16 fr _1 M ; points are connected by 
lines which extend down to the particle mass limit corre- 
sponding to the unsmoothed density field. Squares highlight 
n p (Mf) obtained for a Gaussian filter and are plotted for 
filter masses larger than the minimum halo mass identified 
in each box. 

The quantity n p (M/) can also be calculated analytically 
from the statistical properties of local maxima for a Gaus- 
sian field (BBKS). The comoving number density of peaks 
of arbitrary height is given by 



... . 29-6^6 
MMf) = 53/32(2*)^' 



(1) 



where R* = v / 3(o"i/c"2) is the characteristic comoving 
length-scale of the peaks, and the <7iS are spectral moments 
of the (smoothed) linear power spectrum P(k): 



1 

2^2 



W(k,M f f P(k)k 2(i+1) dk. 



(2) 



The three curves in Figure 1 show the predicted co- 
moving number density of peaks for various choices of filter 
function and upper limit of integration in eq. (2). The dot- 
dashed curve assumes a Gaussian window function and in- 
tegrates the full linear power spectrum used to generate the 
initial conditions of the simulations. The long-dashed curve 
shows the effect of smoothing with a tophat filter and trun- 
cating the integration at the Nyquist frequency of the 1200 
h^ 1 Mpc box In this case, the predicted number density of 
peaks at large filter masses exceeds that measured directly 
from the smoothed density fields by a significant amount. 
For example, at a filtering scale of ~ 10 15 /i _1 Mg one ex- 
pects roughly an order of magnitude more peaks than are 
actually found in the linear density field. 

The origin of this discrepancy is made clear by con- 
sidering the effect of truncating the integral in eq. (2) at 
fcmax = 27r/R/ (shown as a dotted-line in Figure 1). This 
truncation is chosen explicitly in order to ignore the con- 
tribution to n p (Mf) from power on spatial scales below 
the characteristic size of the filter. Although this effec- 
tively changes the tophat filter function by truncating it in 
k— space, we have verified that imposing the same restric- 
tion on the filter prior to smoothing the initial density field 
of the simulations yields peak number densities consistent 
with what is expected from eq. (1). In the remainder of the 
paper we limit our discussion to results obtained with the 
full tophat filter. 



t The truncation is necessary to avoid the divergence of <72 as 
femax — ► for a tophat filter, since for this choice W(k,Mf) 2 oc 
k~ 6 . 
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Figure 2. Smoothed linear overdensity, <5 P , extrapolated to z = 
for peak particles belonging to a single halo plotted as a function 
of smoothing scale. Red points denote 8 p (Mf) for all peak parti- 
cles associated with the object, defined as those that end up in 
the halo by z = 0. The solid (blue) square shows the peak iden- 
tified on the smoothing scale nearest to the actual mass of the 
halo, with the latter shown as a downward pointing arrow. For 
cosmetic purposes we have randomly perturbed the peak masses 
by one smoothing bin width for all peaks identified on scales 
Mf < 1O 13 /i- 1 M0. 



3.3 Associating density peaks with collapsed 
haloes 

Having identified peaks on the sequence of smoothed density 
fields, 5 S , we next associate peak particles with haloes in 
our friends-of-friends catalogs. This is achieved by scanning 
the subset of particles belonging to each halo for matches 
in the list of peaks. By performing the matching at each 
smoothing scale we thereby build lists of peaks associated 
with each object. The outcome of applying this procedure 
to the most massive halo in our 150 h~ Mpc box is shown 
in Figure 2. Here red points show the linearly extrapolated 
peak overdensity, 8 P , as a function of smoothing scale for all 
peak particles associated with the halo. 

Clearly well resolved objects will contain an abundance 
of peak particles with a variety of characteristic masses. In 
the context of peaks theory, the majority of these peaks 
correspond to the progenitors from which the halo formed. 
Although these peak tracers may be of interest for study- 
ing the distribution and dynamics of an early generation 
of accreted (and possibly disrupted) sub-systems (Diemand 
et al., 2005), the vast majority cannot be responsible for the 
collapse of the halo as a whole. If one is to trust the primary 
tenet of the peaks formalism then a single peak, identified 
on the characteristic smoothing scale Mf m Mfof, should 
be singled out. For the example in Figure 2 this peak is 
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shown as a solid blue square, and the actual halo mass as a 
downward pointing arrow. 

The results presented in the following Sections are based 
on the association of a single and unique peak with each 
FOF halo, which we refer to as the "main" halo peak for 
convenience. For each object we define the main peak as the 
one identified on the smoothing scale for which the filter 
mass, M/, is nearest to the true halo mass, that also ends 
up in the halo at z — 0. This method of associating haloes 
with peaks in the smoothed linear density field provides a 
robust lower limit to the fraction of haloes of mass M that 
are associated with peaks in the density field when smoothed 
on same characteristic scale, and thus provides a strong test 
of the primary assumption of the peaks formalism. 

One shortcoming of this method was highlighted by 
Katz et al. (1993), who found that the complex merger histo- 
ries of some haloes may lead to the ejection of peak particles 
during the formation process. Since our runs lack the tem- 
poral resolution to properly track these events, we asses this 
effect in the following way. If at z = a halo contains no 
peak particle from the appropriately filtered linear density 
field (the one with M/ ~ Mfof), we trace the halo parti- 
cles back to the initial redshift and search for one such peak 
within one smoothing radius Rf = (3Mpof/47tp) 1 ^ 3 from 
the proto-halo's center of mass. If a peak particle is found 
to lie within R/ at the initial redshift then we attribute the 
formation of that halo to the peak, and include the particle 
in the peak catalogs. This typically increases the fraction of 
collapsed haloes that are properly identified with peaks in 
the smoothed density field by 5—10%, with a slight tendency 
to be higher for lower mass objects. 



4 RESULTS 

4.1 Peak versus halo masses 

For all haloes containing peak particles we plot the esti- 
mated main peak mass, M pca k, versus the actual halo mass 
in Figure 3. To avoid discreteness effects in M pea k we grid 
the data and use contours to highlight the 99.9, 98 th , and 
90 th percentiles of the distribution. For clarity, the 99.9 and 
98 th percentiles are shown as filled grey contours for our 
150 h' 1 Mpc, but left open for the 1200 h' 1 Mpc-box run. 
Blue curves are used for the 150 hT 1 Mpc box; red curves 
show results for our 1200 h' 1 Mpc box. The solid white 
line corresponds to M pca k = Mfof, and dashed lines show 
the effective resolution of our sequence of smoothed density 
fields. 

Although the scatter in M poa k at fixed halo mass is sig- 
nificant, we find that the majority of haloes can be properly 
associated with peaks of the same characteristic mass in the 
smoothed linear density field. This trend is only weakly de- 
pendent on halo mass, implying that the majority of dark 
matter haloes identified in cosmological simulations do, in 
fact, collapse around density maxima in the linear density 
field. This can be seen in the strong correlation between the 
two masses at the 90 per cent level, which follows very closely 
the 1:1 line anticipated by peaks theory. Finally, it is worth 
noting that the results of our two runs agree remarkably well 
at the 98 th percentiles; even down to the 32— particle halo 
limit of our group finder the two distributions are virtually 
indistinguishable. 
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Figure 3. Main peak mass plotted as a function of the true 
friends-of-friends mass for all haloes in our cosmological simu- 
lations that contain peaks. Peak mass is determined by search- 
ing the halo particle set for peak tracers identified on the lin- 
ear density field smoothed with a top-hat filter contained mass 
ss Mfof • For cases where no such peak particle exists we scan 
the Lagrangian region within the radius Rr = (3Mfof/47tp) 1,/3 
from the proto-haloes center of mass for such a match, or identify 
peak tracers from smoothing scales that most accurately represent 
the true halo mass. For our 150 h^ 1 Mpc box the light shaded re- 
gions enclose the 99.9 percentiles of the distribution; dark shaded 
regions the 98 th percentile, with filled color contours highlighting 
the 90 th percentile. Results from our 1200 h- 1 Mpc box have 
been over-plotted using (unshaded) red contours. 



These results are shown in another way in Figure 4, 
where we plot the cumulative fraction of haloes of mass 
greater than M that also contain peak particles. Different 
line-styles show different subsets of haloes defined in terms 
of the "best" peak mass; solid curves to the fraction of haloes 
for which M pca k ~ Mfof (i.e., those for which the peak mass 
lies within one smoothing bin- width of the true halo mass) , 
and dot-dashed lines to a slightly more relaxed criterion: 
1/2 < MFOF/M pea k < 2. For comparison, we show the cu- 
mulative fraction of haloes associated with peaks from any 
smoothing scale with dashed lines. 

From Figure 4 it is clear that virtually all haloes re- 
solved with more than 100 particles contain at least one 
peak particle. This is in stark contrast with the results of 
Katz et al. (1993), who found that at least 30% of haloes 
in their cosmological simulations to contain no peak tracer 
at all. The large discrepancy between our results and those 
of Katz et al. (1993) is mainly due to the fact that these 
authors matched haloes with peaks identified on only two 
smoothing scales: M/ ~ 2 x IO^Mq and ~ 8 x IO^Mq 
(corresponding to filter volumes containing 9 and 36 parti- 
cles for their simulations). Here we associated haloes with 
peaks on a sequence of smoothing scales spanning the entire 
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Figure 4. Cumulative fraction of haloes above a given mass 
that also contain peak particles. Dashed curves correspond to 
the fraction of haloes that contain peaks from any smoothing 
scale; solid curves with (Poisson) errors to the cumulative frac- 
tion of haloes associated with peaks of the same characteristic 
mass. Dot-dashed curves assume a slightly less restrictive crite- 
rion, showing the fraction of haloes whose "main" peak lies in the 
range 1/2 < MFOF/M pea ]j < 2. As in other Figures, we use blue 
curves to distinguish the results of our 150 h~ 1 Mpc box from 
those of our 1200 h" 1 Mpc box, which are shown in red 



mass range of collapsed objects in our simulations and find 
that Si 99% of haloes with Nfof > 100 particles contain 
at least one peak particle. In the remainder of the paper we 
consider only those haloes that contain at least 100 particles 
unless explicitly stated otherwise. 

Furthermore, the majority of the haloes also contain a 
peak particle from the density field smoothed on the halo 
mass scale, implying that they form in the vicinity of peaks 
of the same characteristic mass. In our 150 ft" 1 Mpc box, 
~71% of haloes form from such a peak. This fraction is 
slightly higher, at ~79% for haloes in our 1200 ft" 1 Mpc 
box. By scanning a range of smoothing scales for peak 
tracers this fraction increases; £ 84% of haloes in our 
150 ft" 1 Mpc box have a peak in the smoothing range 
Mpeak/2 < Mfof < 2 x M pcak , and ~90% of haloes in our 
1200 ft" 1 Mpc box have such a peak. 



4.2 Impact of accretion geometry on estimates of 

Mp ca k 

Although the majority of haloes in our simulations collapse 
near local maxima of similar characteristic mass in the linear 
density field, a significant fraction do not. For example, in 
our 150 ft -1 Mpc box ~16% of proto-haloes contain no peak 
tracer in the mass range Mfof /2 < M pca k < 2 Mfof • These 
are not exclusively low-mass systems; of the subset of haloes 
with > 10 4 particles (or a mass of > 2.43x lO 12 ft _1 M ) 



we find that ~ 5.6% contain no peak in the quoted mass 
range, and a few (~ 1%) can only be matched with peaks 
whose mass differs from the true halo mass by at least a 
factor of 4. Figure 4 shows that similar results hold for 
our 1200 ft" 1 Mpc box run, in which ~ 3.5% of all haloes 
can only be matched with linear density peaks on scales 
M/ Mfof/4. Although these haloes do, in fact, host peak 
tracers we hereafter refer to them as "peakless" haloes for 
convenience (the term reflects the discrepancy between their 
appareent peak mass and the actual halo mass scale). 

The top two panels of Figure 5 show examples of two 
peakless proto-haloes found in our 150 ft" 1 Mpc box. These 
were selected at random from the subset of peakless haloes 
that also have Mfof > lO 13 ft _1 M0 (note that these haloes 
are well resolved and contain ~ 11 x 10 (left panel) and 
~ 5.5 x 10 4 (right panel) particles). The choice of color 
scheme highlights the density contrast field, 1 + 8 S , in the 
halo vicinity after smoothing on the halo mass scale. Or- 
ange curves correspond to a density contrast 5 S = 1, and 
yellow curves to the density contrast associated with the 
spherical collapse barrier, S sc = 1.686. Regions of higher 
density contrast are shown as thin white lines to guide the 
eye. For comparison, the lower two panels plot examples of 
two proto-haloes for which M pea k ~ Mfof- In each panel, 
red dots show the projected positions of the haloes relative 
to the surrounding density field. 

One thing to note about these projections is that the 
characteristic filter radii for the peakless haloes, defined as 
R/ = (3A//47r7j) 1//3 , differ substantially from the character- 
istic sizes of the collapsed regions (proto-haloes have very 
different sizes along their principal axes). For the halo shown 
in the upper left-hand panel, Rf — 4.37 ft" 1 Mpc, yet the 
median and maximum Lagrangian distance to halo particles 
(measured with respect to the proto-halo's center of mass) 
are 4.8 ft" 1 Mpc and 14.8 ft" 1 Mpc, respectively. This sug- 
gests that the characteristic filter radius is not a good esti- 
mate of the size of proto-haloes in the linear density field. For 
the halo in the upper right-hand panel the situation is simi- 
lar; Rf = 3.46 ft" 1 Mpc for this object, which is significantly 
smaller that the true linear dimensions of the collapsed re- 
gion, which spans £ 9 ft" 1 Mpc. This suggests that sim- 
ple spherical filtering masks the complex processes by which 
dark matter haloes form, and may lead to complications in 
associating peaks in the linear density field with collapsed 
objects. 

We explore these difficulties further in Figure 6, which 
examines the distribution of peak particles in and around 
the two peakless proto-haloes shown in the upper panels of 
Figure 5. The top panels show, using red points, the mass 
dependence of peak overdensities for peak tracers belonging 
to each FOF halo. For comparison, open (blue) squares plot, 
as a function of M/, the overdensity of the peak nearest to 
the proto-haloes center of mass as the smoothing scale is 
varied, without requiring that this peak end-up in any halo 
by z = 0. The distance to the nearest peak from the proto- 
haloes center of mass, r scp , is shown in the lower panel. 

Figure 6 exposes the intrinsic difficulties with assign- 
ing peaks to collapsed objects in cosmological simulations. 
Consider first the results in the left-hand panels. Despite 
containing over 1.1 x 10 particles, this halo hosts no peak 
tracers from any filtering scale M f ^ Mfof /4. This is clearly 
at odds with the naive expectations of peaks theory. Fur- 
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Figure 5. Examples of the overdensity field in the vicinity of four haloes found in our 150 /i -1 Mpc box. Top panels show haloes for 
which Mpcak < Mpop/4; bottom panels show two similar mass haloes with Mp,,^ ss MpoF- Linearly extrapolated density fields have 
been smoothed with a tophat filter containing mass MpoF- Particles that make-up the FOF haloes at z = are show as red dots. In all 
panels contours are used to highlight the density gradient in the neighborhood of the halo; a density contrast <5 S = 1 is shown as a orange 
curve; the threshold for spherical collapse, S s = 1.686, is shown as a thick yellow line. Note that these haloes were randomly selected 
from those having Mfqf 10 ls h~ 1 Mg, and all contain more than 5 X 10 4 particles. 



ther scrutiny, however, reveals the presence of a peak of the 
correct characteristic mass in the immediate vicinity of the 
proto-halo which, because of its complex accretion geometry, 
never becomes part of the system. Furthermore, because this 
peak is separated from the proto-halo by more than one filter 
radius it is never associated with the collapse of this object, 



in spite of its proximity to the proto-halo's center of mass. 
The "jump" in the value of 5 P and r sep at M pca k w Mfof 
betrays the fact that, as R/ is increases, isolated peaks may 
merge with others in their vicinity which can be problematic 
when attempting to associate collapsed objects with peaks 
in S s . 
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Figure 6. Top panels show the distribution of peak overdensities, <5 P , as a function of smoothing scale for peak tracers found in the two 
proto-haloes shown in the upper panels of Figure 5. Red points show all peak particles that, by z = 0, are part of the FOF halo; blue 
squares show the peak identified on each 8 S that lies closest to the proto-halo's center of mass, regardless of its association to the halo. 
Downward pointing arrows mark, in each case, the halo mass, MpQF, and the nearest filtering mass scale for which peak tracers are 
found in the halo. Bottom panels show the comoving separation between the proto-halo's center of mass and the nearest peak identified 
on each smoothing scale. For comparison, we show the filter radius Rf(Mj-) as a solid black line; the dot-dashed line shows the mean 
inter-peak separation as a function of smoothing scale. 



The right-hand panels of Figures 5 and Figure 6 
clarify the situation. This object, with a total mass 
of ~ 1.3 x 10 li h~ 1 M(r ) , contains no peak tracers 
from any smoothing scale in the mass range 2.1 x 
10 12 /i _1 M Q <M f < 10 14 /i _1 M Q , corresponding to a differ- 



ence in peak and halo mass of at least a factor of ~ 6. 
However, by tracing the trajectory of the peak closest to 
the proto-halo's center, it is clear that this object does, in 
fact, collapse in the neighborhood of a density maxima of 
the appropriate characteristic mass in the smoothed density 
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field. The upward pointing arrow marks the smoothing scale 
closest to Mpof on which a peak lies within the character- 
istic filter radius of the halo's center of mass. Although this 
peak particle never actually ends up in the halo, it would be 
illusory to ignore it when attempting to map proto-haloes 
to peaks in the initial conditions. 

This emphasizes prior suspicions that simple spherical 
filtering of the linear density field can be misleading, and 
highlights the dangers of over-interpreting the failure of pre- 
vious work to provide a convincing one-to-one mapping be- 
tween haloes and peaks in 8 S . None the less, it is important 
to emphasize before proceeding that not all haloes form in 
the immediate vicinity of similar scale peaks. We show this 
in Figure 7, which plots the distribution of minimum sep- 
arations between proto-halo centers and peaks of the same 
characteristic mass. Different colored curves correspond to 
haloes of different mass, selected to lie in several narrow 
mass bins (the bin width, in each case, coincides with the 
smoothing bin width for our sequence of smoothed density 
fields) . We show only results for haloes in our 1200 h' 1 Mpc 
box. For each mass bin, the percentage of haloes whose 
nearest peak lies within one filter radius is indicated. As 
expected from Figure 4, high mass haloes virtually always 
reside near peaks of the same characteristic mass. Towards 
lower masses, however, the distribution of minimum sepa- 
rations becomes increasingly bimodal; one "hump" corre- 
sponds to the majority of haloes that reside near peaks ini- 
tially, the other is induced by the boundary imposed by the 
mean inter-peak separation. The vertical lines in Figure 7 
show the average separation between peaks on each mass 
scale and clearly show that the vast majority of haloes lie 
within this characteristic distance scale from their nearest 
peak, even if they are physically unassociated with it. 

In the next section, we examine some of the statisti- 
cal properties of peakless haloes, which we contrast directly 
with the bulk of the halo population. To this end, we di- 
vide the sample of haloes into two subsets: those clearly 
associated with peaks in <5 S , for which M pca k « Mfof, with 
which we directly contrast those whose best peak mass is 
M peak < Mfqf/4. 




log r, ep / R f 

Figure 7. Distribution of minimum separations between proto- 
halo centers-of-mass and peaks of the same characteristic mass for 
haloes in our 1200 h~ 1 Mpc box. Different lines styles and colors 
correspond to haloes of different mass, which are selected to lie 
in a narrow mass bin of logarithmic width AlogMpQF ~ 0.129 
(equal to the logarithmic spacing between our set of smoothed 
linear density fields). To aid in the comparison between the var- 
ious haloes samples we have normalized the separations by the 
characteristic filter radius for each mass. Curves are labeled by 
the mass-bin mid-point used to identify each set of haloes, and 
correspond to the filter mass-scale on which the density field was 
smoothed prior to identifing peaks. For each mass scale we also 
indicate the fraction of haloes that lie within one filter radius of 
the nearest peak. The vertical lines, which adopt the same color 
coding, show the mean inter-peak separation for each smoothed 
field. Clearly high mass haloes initially reside near peaks of the 
same mass in the linear density field, whereas for low-mass sys- 
tems this is not always the case. 



4.3 The importance of velocity shear 

Another point worth mentioning is that both haloes shown 
in the upper panels of Figure 5 appear to collapse along 
high-density ridges that connect massive structures in the 
linear density field. This is, in fact, a common occurence for 
peakless haloes found in our simulations, and is not unique 
to the examples plotted in Figure 5. The proto-halo shown 
in the upper-left panel of Figure 5 is initially elongated in 
the direction transverse to the density ridge. This is consis- 
tent with the ellipsoidal collapse model in which maximum 
compression occurs along the primary axis of the shear field, 
which is presumably orthogonal to the ridge. 

A possible implication is that the external tides associ- 
ated with the surrounding large scale structure may have a 
non-negligible impact on the collapse and subsequent virial- 
ization of these systems (e.g. Lilje et al., 1986). As discussed 
by van de Weygaert & Babul (1994) (see also Hoffman, 1986; 
Bertschinger & Jain, 1994), a strong external shear is able 
to accelerate the growth of a small scale fluctuations leading 
to the formation of massive haloes in the absence of a cor- 



responding peak of similar scale in the linear density field. 
This is because convergent velocity flows induced by large 
scale structure may enhance the collapse of nearby small 
scale objects, an effect which may play an important role 
in the formation of low mass haloes (Hoffman, 1986, 1989). 
Tidal forces are also capable of altering the angular momen- 
tum acquisition of proto-haloes (e.g. White, 1984; Porciani 
et al., 2002a), as well as their kinematic and spatial distri- 
bution (e.g. Binney & Silk, 1979; Dubinski, 1992). 

We quantify the linear tidal field at proto-halo centers 
in terms of the deformation tensor 



Di 



G> 2 $ 

dxidx; 



(3) 



Since the trace of Dij is equal to 5, we define the velocity 
shear as its traceless part, Tij = Dij — (Du/3)5ij. Here 5ij 
is the Kronecker delta (not to be confused with the density 
contrast), and $(x) is the peculiar gravitational potential 
which is related to 8 through Poisson's equation, 

V 2 $ = 5. (4) 
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Figure 8. Prolateness and ellipticity of the velocity shear ten- 
sor measured at proto-halo centers in both of our simulations. 
Only haloes with N > 100 particles are shown. Blue contours 
show the distribution for proto-haloes properly matched with 
peaks in the density field when smoothed on the same character- 
istic scale. Green contours show the corresponding distributions, 
but for proto-haloes having Mpcak Mpop/4. The corresponding 
probability distributions are shown in the upper and right-hand 
panels. 



In general, the collapse time depends not only on the 
trace of Dij (i.e., on 5) but rather on its three eigenval- 
ues, which determine the growing mode perturbations: local 
maxima in 5 are not necessarily local minima in collapse 
time. We quantify the tidal field at a given point in terms 
of the eigenvalues of £)y, Ai > A 2 > A3, which determine 
whether the initial flow around a given point is of compres- 
sion or dilation. Equivalent quantities are the shear elliptic- 
ity, e, prolaticity, p, and density contrast, 5. For 5 > 0, we 
have e > and — e < p < e; the relationship to the eigen- 
values of the tidal tensor are conventionally given by (e.g. 
BBKS) 



Ai — A3 



25 



and 



Ai + A3 — 2A2 
25 ' 



(5) 



(6) 



It is customary to quantify the velocity shear at a given 
point in terms of the modified eigenvalues: ti = A,: — 5/3. 
The sign of the U's determine whether local flows promote 
collapse around a given point (positive), or impede it (nega- 
tive). Since ^ . ti = 0, we may also characterize the velocity 
shear in terms of its ellipticity, e s , and prolateness, p 3 , where 



(7) 



h-t 3 3(ti+i 3 ) 
e s = , and p s = , 

and 00 is given by eq. 2. The magnitude of p a determines 



Figure 9. Intermediate-to-major (b/a) versus minor-to- 
intermediate (c/b) axis ratios for the same subsamples of proto- 
haloes plotted in Figure 8. Perfectly oblate haloes have b/a = 1, 
while perfectly prolate objects have c/b = 1; "maximally triaxial" 
haloes have b/a = c/b. Blue contours highlight the regions occu- 
pied by the subsample of haloes having M pca k ss Mpop; green 
contours show the proto-haloes with Mpe^ ^ Mfqf/4. 



the relative contributions from velocity shear along the first 
and third axis of the flow; the sign of p s indicates whether 
the flow along the intermediate axis is inward (p s < 0) or 
outward (p s > 0). 

In Figure 8 we plot the distribution of the ellipticity, e s , 
and prolaticity, p s , of the velocity shear measured at the cen- 
ters of mass of two subsamples of proto-haloes after stacking 
all (Nfof > 100) haloes in both of our simulations. For all 
haloes e s and p s are measured from the tidal field smoothed 
on the mass scale of the halo (note that we have normalized 
by the mass variance, cto(M), to minimize any subtle mass 
dependence in the trends). Blue curves show the results for 
all proto-haloes properly identified with peaks of the same 
mass in the linear density field; green contours to peakless 
haloes. Curves mark the isodensity contours which enclose 
the 98 th , 90 th , 50*\ and 25 th percentiles of the distribution. 
The top and side panels show the corresponding probability 
distributions for e s and p s . 

The velocity shear measured at proto-halo centers is 
slightly but significantly different for the two halo subsam- 
ples. Peakless haloes have slightly more negative prolaticities 
implying that, on average, ti ^ (since ti carries the oppo- 
site sign to p s ). Recalling that . tj = 0, and £3 < 0, this 
implies that the magnitude of the velocity outflow along the 
direction of £3 is, on average, larger than the magnitude of 
the velocity inflow along the primary and intermediate axes 
of the shear field, i.e. \tz \ Si \ti\. This indicates a slightly 
higher tendency for the velocity shear to aid in the collapse 
of these density peaks along two, rather than one axis of the 
shear field. Considering all proto-haloes with Nfof > 100 
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particles, we find that ~52% of those with M pca k ~ Mpop 
have t2 > 0, whereas ti > for ~64% of peakless haloes. 

The differences in velocity shear for the two subsamples 
of haloes highlighted in Figure 8 suggests that the geometry 
of the Lagrangian region from which these haloes grow may 
also differ systematically. For example, a halo with ti > 
t% > is expected to form from a flatter initial configuration 
than a similar mass object with ti > and < 0, since in 
the latter case the velocity shear impedes collapse along the 
intermediate axis rather than aids it. 

We show this explicitly in Figure 9, which plots the 
intermediate-to-major (b/a) versus minor-to-intermediate 
(c/b) axis ratios for the same subsamples of proto-haloes 
used above. (Axis ratios are computed from the eigenval- 
ues of the proto-halo's intertia tensor after weighting each 
particle position by its distance from the center of mass.) 
Note that a perfectly prolate spheroid has c/b = 1, whereas 
an oblate spheroid corresponds to b/a = 1; a maximally 
triaxial object has b/a — c/b. Contours highlight the same 
iso-density regions as those in Figure 8. It is clear that proto- 
haloes with Mp Ca k « Mfof (which accounts for more than 
three-quarters of all haloes) evolve from initial Lagrangian 
regions whose geometries are on average triaxial, with no 
obvious preference for prolate or oblate configurations. How- 
ever, as anticipated from Figure 8, peakless proto-haloes 
clearly deviate from the average, showing a much stronger 
tendency to evolve from initially more flattened configura- 
tions, suggesting that these halo subsets experience very 
different accretion geometries. Since the eigenvalues of the 
tidal field and proto-halo intertia tensor are nearly perfectly 
aligned (Lee & Pen, 2000; Porciani et al., 2002b; Lee et al., 

2009) , we speculate that the initial proto-halo shape is dic- 
tated by the large scale tidal field. We will address this issue 
in forthcoming work (Ludlow & Porciani, 2010, in prep.). 

4.4 Assembly bias 

Sheth & Tormen (2004), Gao et al. (2005) and subsequent 
studies have conclusively demonstrated that the clustering 
strength of dark matter haloes depends on parameters other 
than their mass. Halo properties such as formation time, 
concentration, spin, shape, internal dynamics, and substruc- 
ture content all show correlations with the strength of clus- 
tering (e.g. Wechsler et al., 2006; Croton et al., 2007; Bett 
et al., 2007; Gao & White, 2007; Faltenbacher & White, 

2010) . Somewhat suprisingly, it has been noted that the de- 
pendence of clustering on these parameters does not follow 
intuitively from the intrinsic parameter correlations. For ex- 
ample, older haloes are on average more concentrated than 
younger ones, but the dependence of clustering on forma- 
tion time and concentration exhibit an opposite behaviour. 
This implies that the clustering properties of dark matter 
haloes may depend on a more fundamental parameter, or 
on a number of aspects of the variety of halo formation his- 
tories. 

In Section 4.1 we identified two classes of haloes: those 
clearly associated with similar scale peaks in the linear den- 
sity field, and those with no such peak. We have shown 
that the velocity shear and proto-halo shape for these two 
subsamples are systematically different; peakless haloes are 
more strongly sheared and evolve from more oblate initial 
configurations than the average halo. The possibility that 
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Figure 10. Two-point correlation function for proto-haloes span- 
ning two separate mass ranges: 2.4 X I0 10 h~ 1 Mq < MpQF < 
6.1x10 10 /i- 1 Mq (solid circles), and 6.1 xlO lo /i -1 M < M PO f < 
15.3 X 10 10 h- 1 M Q (open diamonds). These mass ranges include 
haloes having between 100 and 250 particles, and between 250 
and 625 particles, respectively. All correlation functions are com- 
puted at the initial redshift of our 150 h~ 1 Mpc-box simulation 
using proto-halo centers-of-mass. The dotted black curves show 
£(r) for all haloes in each mass bin. The solid blue curves dis- 
tinguish haloes for which M poa k fs MpQF from peakless haloes, 
which are shown using dot-dashed green lines. For comparison, 
the spread in filter radii for haloes in each mass range is shown 
as shaded vertical bands. 



these differences may be influenced by environment moti- 
vates us to explore the statistical properties of their spatial 
distributions by means of the two-point correlation function, 
£(r). 

Figure 10 shows £(r) for proto-halo centers measured at 
the initial redshift of one our simulations. To avoid mixing 
results from haloes of widely different masses we plot here 
only haloes found in our 150 h~ l Mpc box, that lie in two 
narrow mass bins: 2.4 x 10 10 < M FO F/(ft~ 1 M ) < 6.1 x 10 10 
(shown as connected dots in Figure 10), and 6.1 x 10 10 < 
Mfof/(/i _1 Mq) < 15. 3x 10 10 (connected diamonds). These 
mass bins correspond to haloes having between 100 and 250 
particles, and between 250 and 625 particles, respectively. 
The dotted line shows the auto-correlation function of all 
proto-haloes in each mass bin; blue and green curves show 
the effect of limiting the calculation of f(r) to haloes for 
which Mp 0a k ~ Mfof, and those with M pca k < Mfof/4, 
respectively. The separations plotted on the horizontal axis 
are in co-moving units. 

Clearly these halo subsamples exhibit very different ini- 
tial clustering properties. On small scales { S^l h^ 1 Mpc) 
peakless haloes are significantly more clustered than aver- 
age, but the trend is reversed for larger separations. The cor- 
relation function of peakless haloes also has a different slope 
than that of the full halo sample. Between ~ 2—10 h~ l Mpc, 
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Figure 11. Auto-correlation functions for the same halo subsam- 
ple used in Figure 10, but plotted at z = 0. Line-styles and colors 
have the same meaning as before. Note the excess clustering of 
peakless haloes on scales 1 — 2 h~ 1 Mpc. 

for example, the amplitude of the correlation functions can 
differ by over an order of magnitude. Although this effect is 
strongest in the lowest halo mass bin it is still clearly present 
in both, suggesting that it is not particularly sensitive to the 
selected halo mass range. We have noted, however, that a 
small but significant fraction of peakless haloes tend to re- 
side in the immediate vicinity of more massive systems in 
our initial conditions. For example, ~10% of peakless haloes 
in the mass range 2.4 x f0 10 < M FO f/(/i" 1 M q ) < 6.1 x 10 10 
lie within the characteristic radius R= (3M/47rp) 1 ^ 3 of a 
halo with M > 10 13 ' 5 h' 1 Mpc; for those halos with M pcak » 
Mfof this fraction is around 3%. Excluding these halos from 
the calculation of £ (r) results in a significant change in clus- 
tering strength of peakless halos on scales £ 3 h~ x Mpc, 
which then more closely follow the blue curves in Figure 10. 

In Figure 11 we plot the z — two-point auto- 
correlation functions for the same halo subsamples, using 
the same linestyles and colors as in Figure 10. At the small- 
est separations (r ~ 2 h~ x Mpc), the difference in clustering 
strength for the two subsamples differs by a factor of ~ 4. 
Combined with the results on the shapes of proto-haloes 
in Figure 9, this suggests an intimate connection between 
the accretion geometry of a given halo and its proximity 
to neighboring systems (see also, Wang et al., 2007; Hahn 
et al., 2009). Surprisingly, we find that the z = correlation 
function of peakless haloes now exceeds that of haloes with 
Mpeak ~ Mfof at all separations. This is particularly in- 
triguing since, on linear scales, the Eurlerian bias 6_b should 
be related to the Lagrangian bias, 6_L,by&_B = 6_L + l (Mo & 
White, 1996). Further study aimed at addressing the phys- 
ical origin of this behaviour is currently ongoing. 

Although the results presented in Figures 10 and 11 
are limited to haloes in a narrow mass range, and only one 



cosmological simulation, they appear robust to the specifics 
of this selection criterion. The fact that the initial shape of 
proto-haloes, and the strength by which they are sheared, 
is related to their clustering strength may have interesting 
implications for galaxy formation models. 



5 SUMMARY 

We have used two fully cosmological simulations of struc- 
ture growth in the standard ACDM cosmogony to study the 
conditions from which dark matter haloes emerge from the 
linear density field. The primary focus of this work was to 
perform stringent yet simple tests of the central ansatz of 
the peaks formalism, which states that dark matter haloes 
of mass M evolve from peaks in the linear density field when 
the latter is smoothed with a filter of the same characteristic 
mass. 

We identify possible sites for halo collapse by identifying 
local maxima, or peaks, on a sequence of smoothed density 
fields that span the entire mass range of FOF haloes in our 
simulations. For each smoothed density field we build a list 
of peak particles by tagging those nearest to each peak grid 
point, which we then cross-correlate with the halo particles 
in our friends-of-friends catalogs. Our main results can be 
summarized as follows. 

• We define the "main" halo peak as the one identified 
on the smoothed field for which M/ is closest to the true 
halo mass, and use this to "predict" the expected mass 
of the halo. We find that as many as ~70% of all haloes 
(with Nfof > 100) in both of our simulations can be prop- 
erly identified with peaks in the linear density field when 
smoothed on the mass scale of the halo. This fraction de- 
pends systematically on halo mass, with as many as ~91% 
of haloes with M> 5 x 1O 14 /i _1 M0 forming in the vicinity of 
peaks of the expected characteristic mass. By scanning each 
halo for peaks over a broader range of smoothing scales we 
find better agreement: more than ~ 85% of haloes form from 
peaks in the mass range M pca k/2 ^ Mfof ~ 2 x M poa k. 

• Although the majority of haloes in our simulations 
form preferentially around similar-scale peaks in the lin- 
ear density field a small but significant fraction show a 
considerable disparity between the predicted and measured 
masses. For example, ~20% of haloes with > 100 parti- 
cles have Mp 0a k < Mfof, of which ~ 15% contain no peak 
tracers from any scale within a factor of four of the true 
halo mass. These are typically low-mass objects, but several 
have M £ lO 13 ft _1 M0 and contain more than 10 4 particles. 
Haloes for which M pca k Mfof/4 we have referred to as 
"peakless" haloes for convenience. 

• By contrasting directly the properties of peakless proto- 
haloes with those of the population bulk we find slight but 
systematic differences in both the linear tidal shear and ini- 
tial shapes of the two samples. The tidal shear (measured 
at proto-halo centers, and smoothed on the halo mass scale) 
acting on the proto-haloes implies that, on average, peakless 
haloes are tidally compressed along two axes of the initial 
tidal tensor, and expanding along one. On the other hand, 
those with M pea k « Mfof are characterised by roughly 
equal compression and expansion along two axes of the tidal 
field, and a null flow along the third. 
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• Peakless haloes also evolve from more oblate initial con- 
figurations than the nearly triaxial geometries characterising 
the mean. Haloes for which M pca k ~ Mpop typically collapse 
from Lagrangian regions which are initially triaxial with typ- 
ical axis ratios b/a ~ c/b « 0.75. Peakless haloes, on the 
other hand, are initially more flattened, having b/a « 0.75 
and c/b « 0.5. The differences in the initial shapes of proto- 
haloes are qualitatively consistent with the differences in 
their tidal shear; we plan to address the interrelationship 
between tidal forces and accretion geometry in a forthcom- 
ing paper (Ludlow & Porciani, 2010, in prep.). 

• The auto-correlation functions of normal proto-haloes 
and of peakless haloes show significant differences at both 
the initial redshift of the simulation as well as at z = 0; on 
scales 1 h~ Mpc peakless haloes are significantly more 
clustered both initially as well as today. The evolution of 
the correlation functions for the different subsamples, how- 
ever, appears quite different; on scales larger than a few 
Megaparsecs peakless haloes are today more clustered than 
average, whereas they are initially less clustered at the same 
separations. Further study is required to address the origin 
of this behaviour. Combined with the differences in the ini- 
tial shapes and shear properties of these samples, these re- 
sults imply that the accretion geometry (and hence the mass 
acquisition history) of highly clustered haloes can differ sub- 
stantially from that of field haloes. This may have important 
implications for understanding the physics of galaxy forma- 
tion in dense environments and for understanding the role 
environment plays in the morphological transformation of 
galaxies. 

The main goal of this work has been to provide a simple 
and comprehensive test of the central ansatz of the peaks 
formalism. In doing so we have demonstrated that the ma- 
jority of dark matter haloes identified in our cosmological 
simulations form in the vicinity of linear density peaks of the 
expected characteristic mass in the simulations initial con- 
ditions, yet a significant fraction do not. Our work sparks 
several interesting questions for future study: what is the 
physical origin of the peakless haloes? What implications 
do these objects have for analytic models of structure for- 
mation rooted in peaks theory? What is the nature of their 
clustering properties and their imprint on the assembly bias? 
We plan to address these and other pertinent issues in forth- 
coming work. 
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